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SUMMARY 


In this note, a Chebyshev matrix collocation method is outlined for the solution 
of the Orr-Sommerfeld equation for the Blausius boundary layer. User information 
is provided for FORTRAN program ORRBL which solves the equation by the OR 
method. 

1. PROBLEM FORMULATION 

In this report, we summarize a numerical procedure that was implemented 
to develop the FORTRAN program ORRBL which solves the Orr-Sommerfeld equation 
for two- and three-dimensional small amplitude disturbances. The program is written 
for the Blausius boundary layer but the method can be easily adapted to other parallel 
shear flows. The current version of the program solves for the temporal eigenvalues 
and a new version for the solution of the spatial eigenvalues is under development. 

According to the linear stability theory, the equations of continuity and 


momentum for parallel shear flows can be written as (Ref. 1) 

i(aui -f l3u 3 ) -r D 112 = 0, (1) 

[ D 2 — (a 2 + /3 2 ) — iaR(u — c)]fli = R(u')u 2 + iaRp, (2) 

[D 2 — (a 2 + (3 2 ) — iaR(u — c)]u 2 = RDp, (3) 

[D~ — (a 2 + (3 2 ) — iaR(u — c)]« 3 = if3Rp. (4) 

Note that these equations satisfy the disturbance field defined as 

v-i = v-i(y) exp[i(ax + (3z) — iact). (5) 


In eqs. (1-5), D = d/dy , u is the mean (Blausius) velocity profile, u is its first 
derivative y, a is the wave number along the streamwise (a:i) direction, /? is the wave 
number along the spanwise (x 3 ) direction. Also, the complex dimensionless amplitude 
functions uj, u 2 , “3 and p are functions of the cross-stream variable ® 2 . For temporally 
evolving disturbances, the eigenvalue c is complex; its real part corresponds to the 



phase speed of the wave and its imaginary part indicates amplification in time, i. e. for 
given Reynolds number (J?) and wave numbers (a and /3) if c/ < 0 for the least-damped 
eigenvalue, then the flow is stable. However if cj > 0, then the flow is unstable. For 
the Blausius boundary layer, R is defined in terms of the displacement thickness and 
free-stream velocity. 

The three-dimensional disturbance equations can be reduced to the Orr- 
Sommerfeld equation by the Squire transformation for which we define (Ref. 1) 

k 2 = a 2 + f3 2 , 

v = u 2 , (6) 

au i + (3ii3 = kit. 

With this substitution, the Orr-Sommerfeld equation reads 

(D~ — a 2 ) v = ia — c)(D 2 — a 2 )v — (a")-?]. (7) 

For the boundary layer, this equation is solved with the boundary conditions. 

v = Dv = 0 at y = 0 and y — » oo. (8) 

In program ORRBL, eqs. (7,8) are solved by a Cliebyshev matrix colloca- 
tion method which results in a matrix equation 


(A — cB)y = 0. (9) 

Here, A and B are functions of the Chebyshev collocation matrix. This matrix is 
explicitly given in Ref. 3. The matrix eigenvalue problem (eq. 9) is then solved by 
the QR method using the IMSL subroutine EIGZC. The mean velocity profile u is 
obtained by solving the Blausius equation which reads 

/'" + \ff" = 0, (10) 
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by a point-by-point fifth and sixth-order Runge-Kutta-Verner method. 

The current version of the program provides very high accuracy with about 
100 Cliebyshev collocation points; in the program this number is specified by the user. 
Note that Cliebyshev polynomials are defined in the interval — 1 < y < 1, so that for 
application to the boundary layer problem, a coordinate transformation is required. 
For this purpose we use (Ref. 2) 


V = 



(ID 


where Y is a user defined parameter which should be 15 < Y <25 for maximum 
accuracy. Note that this transformation allows the inclusion of the boundary points 
and satisfies boundary conditions exactly. 

2. SOLUTION METHOD 

Chebysliev spectral methods have recently been very successful in boundary 
layer stability and transition investigations. Spectral methods are highly accurate 
and require much less discretization points than finite-difference and finite-element, 
methods. The difficulty of their application to the semi-infinite integration domain of a 
boundary layer has been overcome with suitable mapping or transformation techniques. 
The matrix collocation method can be relatively easily implimented into existing codes. 
A brief overview of this method follows. 

The semi-infinite integration domain 0 < y < oo can be transformed into 
the interval 0 < r/ < 1 using the mapping formula (Ref. 2) 


V = exp(-y), (12) 

with the free transformation parameter Y. In the transformed variable 77 , any function 
u can be expanded into a Cliebyshev series 

N 

u (v) = X {lk ' rfc ( T ?)> ( 13 ) 

k = 0 
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with Chebyshev coefficients Uk (called spectrum of u) and the Chebyshev polnomials 


Tk(rj) = cos(kcos 1 77), (14) 

defined in the interval — 1 < 77 < 1. Defining the ’collocation points’ 

Vj = c °s( 7 jf-) ; J = 0..JV, (15) 

the series (13) can be evaluated. With the spectrum Uk known, derivatives of u can 
then be evaluated. 

The application of this method requires the use of the Chebyshev matrix 
D^ p K The coefficients of this matrix are known. This matrix relates the values of 
u at the collocation points T]j to the values of the p-tli derivative of u at these points: 

^ = (16) 

' fc = 0 

We now define 

u(rjj) = 0 ; j > N/2 + 1 = Nkaif, (IT) 

( N ' even) and apply the transformation (12) to the derivative relation (16) at the 
transformed collocation points 

Hi = —Yln(rjj ) ; j = 0...iV^/ - 1. (18) 

Only one quarter elements of the original collocation matrix D have to be used. The 
boundary condition u = 0 at y — ► 00 is satisfied automatically because of the trans- 
formation (17). 
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SAMPLE CASE 

INPUT FILE (INPUT PARAMETERS) 

The code requires an input data file named ORDATA. The assigned values 
for the following parameters should be included in this file. 

R - (REAL) Reynolds number based on displacement thickness and the free- 
stream velocity. 

a - (REAL) Streamwise wave number. 

f3 - (REAL) Spanwise wave number. For two-dimensional calculations /3 should 
be equal to zero ( f3 = 0). 

Y - (REAL) This parameter is used for the coordinate transformation. A high 
value of Y means mesh clustering towards the boundaries. For best results: 
15 < Y < 25 (Y=20 is suggested). 

N - (INTEGER) Number of the collocation points in the whole domain. It 
should be an even number which is < 100. Due to the coordinate transfor- 
mation, only half of the given N is used for the eigenvalue calculations in 
the code. 

IOPT - (INTEGER) Option parameter. 

IOPT=0 - only the eigenvalues are calculated, 
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I0PT=1 - both the eigenvalues and the eigenvectors are calculated, but 
only the eigenvector corresponding to the least damped eigen- 
value is printed. 

EXAMPLE: 

The following parameters are given, 

R = 998, a = 0.308021, /? = 0, Y = 20, N = 100, IOPT = 1. 

For this case, ORDATA has the following form, 

998. 

0.308021 

0. 

20 . 

100 

1 

Note that the input file starts at the first column of the first row. 

PROGRAM OUTPUT 

ORRBL creates two output files named OREIGE and ORVEL, respectively. 
Output file 1: OREIGE 

This file contains the following information, 

- number of the collocation points used in the half domain, 

- Reynolds number 

- Streamwise wave number (a), 

- Spanwise wave number ( f3 ), 

- Y, 

- Base velocity profile (Blausius solution), 
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column 1 - physical y-coordinates, 
column 2 - mean velocity profile, 

column 3 - first derivative of the mean velocity profile, 
column 4 - second derivative of the mean velocity profile, 

- Index of the unconverged eigenvalue (for a detailed explanation, please refer 
to the parameter listing of EIGZC); if this number is zero, all the eigenvalues 
are converged. 

- EIGZC error number (for a detailed explanation, please refer to the param- 
eter fisting of EIGZC); if this number is zero, there are no errors. 

- a list of the eigenvalues (Note that if one of the diagonals of B matrix is 
zero, the corresponding eigenvalue goes to infinity. Therefore, in such cases, 
the real part of the eigenvalue is set equal to 1 x 10 30 and the imaginary- 
part is zero.), 

column 1 - index, 
column 2 - real part, 
column 3 - imaginary part, 

- Index of the least damped eigenvalue, 

- Least damped eigenvalue (CR : real part, Cl : imaginary part), 

- Eigenvector corresponding to the least damped eigenvalue, (printed for y < 
30) 

column 1 - physical y-coordinates, 

column 2 - real part of the eigenvector, 

column 3 - imaginary part of the eigenvector, 

column 4 - (if IOPT=l) phase angle of the eigenvector. 

Output file 2: ORVEL 

This file is used if IOPT=l and contains the following information, 

- Real and imaginary parts of the v perturbation velocity, (printed for y < 30) 
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column 1 - physical y-coordinates, 
column 2 - real part of the velocity, 
column 3 - imaginary part of the velocity, 

- Real and imaginary parts of the u perturbation velocity, (printed for y < 
30) If the calculation is two-dimensional, this velocity is the same as uj in 
eq. (6). Otherwise it gives i in the same equation, 
column 1 - physical y-coordinates, 
column 2 - real part of the velocity, 
column 3 - imaginary part of the velocity, 

A sample output is given on the following pages. 
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OUTPUT FILE 1 : OREIGE 

NUMBER OF COLLOCATION POINTS IN THE HALF DOMAIN- 51 
REYNOLDS NUMBER- .9980000E+03 
STREAMWISE WAVE NUMBER- .3080210E+00 
SPANWISE WAVE NUMBER- .OOOOOOOE+OO 
Y- . 2000000E+02 


BASE VELOCITY PROFILE (BLAUSIUS SOLUTION): 


Y-COORDINATE 


.0000000E+00 
. 9871228E-02 
. 3950442E-01 
. 8895825E-01 
. 1583310E+00 
. 247761 5E+00 
. 3574299E+00 
. 4875596E+00 
. 6384187E+00 
. 8103224E+00 
. 1003636E+01 
. 1218777E+01 
. 1456221E+01 
. 1716504E+01 
. 2000229E+01 
. 2308071E+01 
•2640783E+01 
. 2999209E+01 
. 3 384286E+01 
. 3797063E+01 
. 4238707E+01 
. 4710523E+01 
. 5213969E+01 
. 5750680E+01 
. 6322492E+01 
. 6931472E+01 
•7579956E+01 
. 8270595E+01 
. 9006405E+01 
. 9790840E+01 
. 1062787E+02 
. 1152210E+02 
. 1247889E+02 
. 1350452E+02 
. 1460645E+02 
. 1579358E+02 
. 1707668E+02 
. 1846893E+02 
. 1998668E+02 
. 2165057E+02 
. 2348718E+02 
. 2553151E+02 
. 2783097E+02 
. 3045207E+02 
. 3349219E+02 
. 3710236E+02 
.4153558E+02 
. 4726618E+02 
. 5535902E+02 
. 6921210E+02 
. 6570835E+03 


U 

.0000000E+00 
. 3277814E-02 
. 1311773E-01 
. 2953910E-01 
. 5257354E-01 
. 8226237E-01 
. 1186498E+00 
. 1617681E+00 
. 2116109E+00 
. 2680876E+00 
. 3309542E+00 
. 3997140E+00 
. 4734899E+00 
. 5508795E+00 
. 6298269E+00 
. 7075738E+00 
. 7807826E+00 
. 8459167E+00 
. 8998962E+00 
. 9408823E+00 
. 96886 58E+00 
. 9856925E+00 
. 9944021E+00 
. 9981862E+00 
. 9995284E+00 
. 9999051E+00 
. 9999858E+00 
. 9999985E+00 
. 9999999E+00 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1000000E+01 
. 1 OOOOOOE+Ol 


D1U 


. 3320573E+00 
. 3320573E+00 
. 3320568E+00 
. 3320509E+00 
. 3320209E+00 
. 3319176E+00 
. 3316380E+00 
. 3309943E+00 
. 3296759E+00 
. 3272077E+00 
. 3229096E+00 
. 3158717E+00 
. 3049715E+00 
. 2889696E+00 
. 2667317E+00 
. 2376064E+00 
. 2019234E+00 
. 1614495E+00 
. 1194987E+00 
. 8037511E-01 
. 4811798E-01 
. 2505886E-01 
. 1106876E-01 
. 4032134E-02 
. 1173826E-02 
. 2635151E-03 
. 4379211E-04 
. 5139409E-05 
. 4033041E-06 
. 1985682E-07 
. 5665424E-09 
.0000000E+00 
.0000000E+00 
.0000000E+00 
.OOOOOOOE+OO 
.0000000E+00 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
. OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 


D2U 


.OOOOOOOE+OO 
-. 2686016 E- 05 
- . 4301865E-04 
-. 2181370E-03 
- . 6909488E-03 
-. 1691347E-02 
- . 3516774E-02 
- . 6 529649E-02 
- . Ill 4 6 51E-01 
- . 1780960E-01 
-.2692628E-01 
-. 3875721E-01 
-. 5323555E-01 
- . 6971703E-01 
- . 8671029E-01 
- . 1017255E+00 
- . 1 11 4798E+00 
-. 1127030E+00 
- . 1035333E+00 
- . 8492700E-01 
- . 6100169E-01 
-. 3755105E-01 
- . 1934685E-01 
- . 8126053E-02 
- . 2700904E-02 
- . 6865512E-03 
- . 1282928E-03 
- . 168 3107E-04 
- . 1469160E-05 
-. 8012280E-07 
- . 2523121E-08 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
.OOOOOOOE+OO 
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INDEX OF THE UNCONVERGED EIGENVALUE- 0 
IMSL EIGENVALUE SUBROUTINE ERROR NUMBER- 0 


INDEX EIGENVALUES 



REAL PART 

IMAGINARY PART 

1 

. 5483498E+15 

- . 3065511E+14 

2 

5434929E+15 

4661634E+14 

3 

- . 5849787E+11 

. 9827289E+11 

4 

-.4918929E+01 

- . 1642465E+02 

5 

. 4966318E+01 

- . 16 3 926 3E+02 

6 

8798325E-01 

4258754E+01 

7 

. 2771130E+00 

- . 399 4019Ef 01 

8 

. 1361247E+00 

-.1875335E+01 

9 

. 2862971E+00 

- . 1592164E+01 

10 

. 2887507E+00 

- . 9702207E+00 

11 

.4100225E+00 

- . 7881547E+00 

12 

. 2901083E+00 

- . 2766870E+00 

13 

. 4263477E+00 

- . 5073704E+00 

14 

. 3641256E+00 

. 7959628E-02 

15 

. 5236208E+00 

- . 3776332E+00 

16 

. 4837004E+00 

- . 1922256E+00 

17 

. 5722661E+00 

3031564E+00 

18 

.67 2 6 531E+0 0 

- . 2586152E+00 

19 

. 7609321E+00 

- .2370853E+00 

20 

. 8816834E+00 

- . 24 576 4 3E+00 

21 

. 8449862E+00 

- . 2181525E+00 

22 

. 9081744E+0 0 

- . 1884448E+00 

23 

. 9363043E+00 

- . 1638718E+00 

24 

. 9557013E+00 

- . 1 381578E+00 

25 

. 9690912E+00 

- . 11 55529E+00 

26 

.9784337E+00 

- . 9610866E-01 

27 

. 9848149E+00 

- . 7975591E-01 

28 

. 9892301E+00 

- . 6609220E-01 

29 

. 9922813E+00 

- . 5473425E-01 

30 

. 9944356E+00 

- . 4 5289 12E-01 

31 

. 9959518E+00 

- . 37 43721E-01 

32 

. 9970472E+00 

3089302 E- 01 

33 

. 9978315E+00 

- . 2543808E-01 

34 

. 9984093E+00 

- . 2088029E-01 

35 

. 9988281E+00 

- . 1707562E-01 

36 

. 9991410E+00 

- . 1389604E-01 

37 

. 9993691E+00 

- . 1124572E-01 

38 

. 9995412E+00 

- . 9038025E-02 

39 

. 9996666E+00 

7207606E-02 

40 

. 9997615E+00 

5694493E-02 

41 

. 9998301E+00 

- . 4452860E-02 

42 

. 9998820E+00 

- . 3440249E-02 

43 

. 9999188E+00 

-.2623668E-02 

44 

. 9999463E+00 

- . 1972373E-02 

45 

. 9999652E+00 

1461875E-02 

46 

•9999788E+00 

- . 1 06 9 2 09E-02 

4 7 

. 9999877E+00 

7753646E-03 

48 

. 9999937E+00 

- . 5630750E-03 

49 

. 9999994E+00 

- . 3360325E-03 

50 

.9999973E+00 

4193824E-03 

51 

. 1000000E+31 

. 0000000E+00 
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INDEX OF THE LEAST DAMPED EIGENVALUE= 14 
LEAST DAMPED EIGENVALUES CR- .3641256E+00 Cl- .7959628E-02 


EIGENVECTOR CORRESPONDING TO THE LEAST DAMPED EIGENVALUE 


Y-COORDINATE 


.0000000 

.0098712 

.0395044 

.0889583 

.1583310 

.2477615 

.3574299 

.4875596 

.6384187 

.8103224 

1.0036358 

1.2187771 

1.4562212 

1.7165043 

2.0002291 

2.3080706 

2.6407832 

2.9992088 

3.3842865 

3.7970632 

4.2387071 

4.7105227 

5.2139687 

5.7506798 

6.3224917 

6.9314718 

7.5799564 

8.2705948 

9.0064048 

9.7908396 

10.6278723 

11.5221016 

12.4788863 

13.5045180 

14.6064469 

15.7935801 

17.0766832 

18.4689309 

19.9866788 

21.6505713 

23.4871801 

25.5315075 

27.8309718 


REAL PART 


. 5319612E-12 
. 8757955E-04 
. 1398679E-02 
. 6994729 E- 02 
. 2141506E-01 
. 4933132E-01 
.9 379128 E- 01 
. 1551793E+00 
. 2316440E+00 
. 3204536E+00 
. 4188825E+00 
. 5237591E+00 
. 6303692E+00 
.7325072E+00 
. 8241697E+00 
. 9006787E+00 
. 9579111E+00 
. 9918822E+00 
. 1000000E+01 
. 9824159E+00 
. 9424044E+00 
. 8856714E+00 
. 8188041E+0Q 
. 7476002E+00 
. 6761006E+00 
. 6065899E+00 
. 5401713E+00 
. 4773642E+00 
. 4184559E+00 
. 3636 377E+00 
. 3130394E+00 
. 2667367E+00 
. 2247520E+00 
. 1870563E+00 
. 1535712E+00 
. 1241720E+00 
. 9869082E-01 
. 7692105E-01 
. 5862157E-01 
. 4352197E-01 
. 3132802E-01 
. 2172749E-01 
. 1439634E-01 


IMAGINARY PART 


. 2445888E-12 
- . 8413862E-04 
- . 1237157E-02 
-. 5425582E-02 
- . 1 395592E-01 
- . 2596601E-01 
-. 3830948E-01 
-. 4700520E-01 
- . 4921888E-01 
-. 4451928E-01 
3495078E-01 
- . 2410007E-01 
-. 1533483E-01 
- . 9909034E-02 
- . 6720111E-02 
-. 4265912E-02 
- . 2188051E-02 
- . 7 423 556E-03 
. 0000000E+00 
. 1766201E-03 
. 9336198E-05 
- . 2698427E-03 
4959850E-03 
-.6074399E-03 
-. 6216548E-03 
- . 5829454E-03 
- . 5255261E-03 
-. 4655761E-03 
- . 4082624E-03 
- . 3547914E-03 
-. 3054238E-03 
- . 2602483E-03 
- . 2192842E-03 
-. 1825061E-03 
- . 1498350E-03 
- . 1211515E-03 
- . 9628974E-04 
- . 7 50 5007E-04 
-. 5719527E-04 
- . 4246346E-04 
- . 3056570E-04 
- . 2119917E-04 
- . 1404593E-04 


PHASE ANGLE 


. 4309628E+00 
- . 7653626E+00 
- . 7241957E+00 
- . 6597268E+00 
- . 5775605E+00 
- . 4845121E+00 
3877735E+00 
- . 2941235E+00 
- . 2093628E+00 
-. 1380423E+00 
- . 8324532E-01 
- . 4598122E-01 
-. 2432195E-01 
- . 13 52673E-01 
- . 8153615E-02 
-. 4736295E-02 
- . 2284186E-02 
- . 7484311E-03 
.0000000E+00 
. 1797814E-03 
. 9906785E-05 
- . 3046758E-03 
-.6057431E-03 
- . 8125195E-03 
- . 9194706E-03 
- . 9610203E-03 
- . 9728876E-03 
- . 97-53 0 57E-0 3 
- . 9756397E-03 
- . 9756725E-03 
- . 9756717E-03 
- . 97 567 4 4E-0 3 
- . 9756715E-03 
- . 9756747E-03 
- . 9756711E-03 
- . 9756753E-03 
- . 9756704E-03 
- . 9756763E-03 
- . 9756690E-0 3 
- . 9756785E-03 
- . 9756661E-03 
- . 9756838E-03 
- . 9756592E-03 
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OUTPUT FILE 2 : ORVEL 


REAL AND IMAGINARY PARTS OF THE V PERTURBATION VELOCITY 


Y-COORDINATE REAL PART 


IMAGINARY PART 


.0000000 
.0098712 
.0395044 
.0889583 
. i583310 
.2477615 
.3574299 
.4875596 
.6384187 
.8103224 
1.0036358 
1.2187771 
1.4562212 
1.7165043 
2.0002291 
2.3080706 
2.6407832 
2.9992088 
3.3842865 
3.7970632 
4.2387071 
4.7105227 
5.2139687 
5.7506798 
6.3224917 
6.9314718 
7.5799564 
8.2705948 
9.0064048 
9.7908396 
10.6278723 
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